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Abstract 



We develop a quasiclassical method based on the path integral formalism, 
to study the influence of disorder on magnetooscillations of the density of 
states and conductivity. The treatment is appropriate for electron systems 
in the presence of a random potential with large correlation length or a ran- 
dom magnetic field, which are characterisitic features of various 2D electronic 
systems presently studied in experiment. In particular, we study the system 
of composite fermions in the fractional quantum Hall effect device, which 
are coupled to the Chern-Simons field and subject to a long-range random 
potential. 
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I. INTRODUCTION. 



In this article we develop a quasiclassical method for studying the influence of disorder 
on magnetooscillations in 2D electronic systems. Our original motivation was the problem of 
quantum particles in random magnetic helds which has attracted considerable interest during 
the last few years. Only quite recently realizations of a random magnetic field acting on a 2D 
electron gas have been prepared and investigated. One possible way of generating a random 
magnetic field is to use a type II disordered superconductor with randomly pinned flux lines 
in an external magnetic field as the substrate for the 2D electron gas |I|. Alternatively, one 
may use a magnetically active substrate such as a demagnetized ferromagnet with randomly 
oriented magnetic domains [Q. 

The main interest in the random magnetic field problem derives, however, from effec- 
tive field theories of interacting electron systems, for which the interaction may be shown 
to be mediated by fictitious gauge helds. The hrst example of this class is the gauge held 
theory of high-Tc superconductivity compounds ^-0. There the gauge helds arise as a tool 
to implement the projection of the Hilbert space onto a subspace of states without double 
occupancy of lattice sites. These gauge helds are long-range correlated, with correlation 
function diverging as 1/q^ for wavevector g — >■ 0. This corresponds to a short-range cor- 
related random magnetic held. The second example is provided by the composite fermion 
picture of the fractional quantum Hall ehect USHTlj. There the electrons are represented by 
fermions with an even number of magnetic hux tubes attached. At half-hlling of the lowest 
Landau band, the average held due to the hux tubes cancels the external magnetic held, 
leaving the problem of fermions moving in a random magnetic held generated by the hux 
tubes. The correlations of the corresponding vector potential are also long-ranged. 

For a short-range correlated magnetic held, implying a long-range correlated vector po- 
tential, the total scattering rate of a charged quantum particle diverges in Born approxima- 
tion P,|^. This is due to the divergence of the diherential scattering cross-section for forward 
scattering. The usual self-consistent treatment of strong scattering has been shown to be 
insufficient 0, as it violates gauge invariance. The dominance of small angle scattering 
calls for a quasiclassical description. We therefore employ the quasiclassical approximation 
method for the path integral representation of this problem. An additional advantage of the 
path integral formalism is the explicit gauge invariance, in contrast to the usual perturbation 
theory methods. 

For the 2D GaAs — AlGaAs heterostructures on which the fractional quantum Hall 
effect (FQHE) is observed, the quasiclassical treatment is appropriate for the following 
resasons. In these systems, the donors are located in a remote layer separated by a large 
distance ~ 50 -i- 80nm from the electron gas plane. Thus, the random potential created 
by these impurities has a large correlation length ~ d^, and the small-angle scattering 
dominates, which can be properly described within the quasiclassical approximation. In a 
strong magnetic field, such that the Landau level filling factor z/ is close to 1/2, a statistical 
transformation can be applied as mentioned above, converting the electrons into so-called 
composite fermions p|-jl0|. Then a random magnetic field appears, on top of the smooth 
random potential. The quasiclassical path integral approach allows to study the effect of 
both types of random field on equal footing. 

The outline of the paper is as follows. In section 2 we use the path integral formalism 
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to calculate the total and the transport scattering rates in long-range correlated random 
potentials. We check that the obtained results are in agreement with the conventional 
perturbation theory. In section 3 we study the influence of a long-range random potential 
on the oscillations of the density of states and the conductivity in a magnetic held. In section 
4 we generalize the results of the two preceding sections on the case of a random magnetic 
held (rather than random potential). Finally, in section 5 we apply the above formalism to 
the system of composite fermions in the FQHE device near v = 1/2, where both long-range 
random potential and random magnetic held are important. Section 6 contains a discussion 
of the results and conclusions. 

Some of the results of this article have been published in the form of short communica- 
tions |1T^,|I5 



II. SCATTERING RATES IN LONG-RANGE RANDOM POTENTIAL FROM 

THE PATH INTEGRAL FORMALISM. 



In this section we develop the quasiclassical path integral formalism and apply it to 
the calculation of total and transport scattering rates in a random potential with large 
correlation length. These results can also be obtained within conventional perturbation 
theory, so that this section has mainly methodological character. It is instructive, however, 
to see how the quasiclassical treatment reproduces results of the perturbation theory, in the 
case when both are applicable. 

We consider a quantum particle of mass m and energy E moving in 2D in a static random 
potential U{r) with Gaussian distribution characterized by the correlation function 

{U{r)U{r')) = W{\r-r'\) (1) 

We assume hF(|r|) to be a smooth function of r depending on the absolute value r = \r\ 
only, with a characteristic length scale ^ 3> A, where A = is the wave length 

(we set h = 1 throughout the paper). We review first results of the standard perturbation 
theory. The total scattering rate can be found in the Born approximation as 

1/ts = ‘2t^ J {dpi)W {\p - p-i_\)5{ep - ep^) , (2) 



where {dp) = d^p/(27r)^; = p^ /2m and 



W{p) = / hF(r) exp(— ipr) 



Eq.(^ can be transformed as follows 



1/ts = 2ttN{E) j ^W{2pFsin^) 
JO 27T 2 



= 2nN{E) J d</) J drrJo ^ppr sin ^ W{r) 

fOO 

= {2n)^N{E) / drrJ^{pFr)W{r) , 

JO 



( 3 ) 



( 4 ) 
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where N{E) = m/27r is the density of states (DOS) on the Fermi snrface and pp = {2mEy^‘^ 
is the Fermi momentnm. The above assnmption of the long-range character of PF(r), ^ 3> A, 
allows to approximate the Bessel fnnction in (^) by its asymptotic expression, yielding 
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1/ts = — / drW{r) , 
Vp Jo 



( 5 ) 



with Vp = pp/m being the Fermi velocity. This resnlt holds nnder the condition of applica- 
bility of the Born approximation, which is 



VfTs > ^ 



( 6 ) 



Expression for the transport scattering rate 1/rtr differs from that for 1/ts by an addi- 
tional factor (1 — COS0). We have thns, in fnll analogy with eq.(|). 



r2n 



d(j) 



1/Ttr = 27 tN{E) / —W{2ppsin-){l-cos(j)) 
Jo Zn / 



= 27tN{E) j d(j) J drr Jf) ^ppT sin ^ {1 — cos (j))W{r) 

poo 

= {2-KfN{E) / drrUhppr) - Jl{ppr)]W{r) . (7) 

Jo 

Using now the relations Jq{z) = —Ji{z)] [zJi{z)]' = zJq{z) and integrating by part, we 
rednce eq.(|^) to the form 



l/^ir = 



27T" 

pI 



N{E) [ dr£{ppr)^r^W{r). 
J dr dr 



( 8 ) 



For a long range potential, 3> A, we hnd thns 






m 



poo 


\ ^ IU'(r)l 


/ dr 


W"(r) + ^ ^ 


Jo 


r 



m 
Pp Jo 



dr 



W\r) 



( 9 ) 



In the last line we nsed the assnmption that W{r) is an analytic fnnction of r at r = 0, so 
that hF'(O) = 0. 

In the remaining part of this section we rederive eqs.(^, (|^) in the path integral approach 
I^ . To calcnlate r^, we consider the (retarded) single-particle Green fnnction 



GR{yR-T) = 



rr{T)=R 

/r(0)=0 



Dr{t) exp lij 



dt 



mr'^ 



U{r) 



( 10 ) 



After averaging over the disorder (which we will denote by angnlar brackets), we get 

I fT rT 



rr{T)=R I r 

{Gr{0,R]T)) = / Dr{t)expli / 

./r(0)=0 1 Jo 






( 11 ) 
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In the absence of the second term in the exponent, eq.(^) describes the free particle Green 
function (G^^(0, R] T)), which can be found exactly by the saddle-point method. The saddle 
point is given by the classical trajectory r{t) = vt] v = R/T, which yields in 2D 

= ( 12 ) 

We use now the fact that the correlator W (r) is smooth, and thus the full expression ( pT|) 
can still be evaluated quaiclassically. Furthermore, we assume the random potential to 
be relatively weak, so that the second term in the action in (|^) can be considered as a 
perturbation. Formally, this means Ets S> 1. Then the effect of the random potential term 
is given by its value on the saddle-point trajectory. Assuming now vpTs S> as in eq.(^), 
we get 



{GUR-,T)) = (Gg’(0,fi;r))exp (-T f°° —W(r)] 

Y Jo Vp J 

= (Gg)(0,i?;T))exp(-T/2r,) , (13) 

with Ts as in eq.(|). We hnd therefore that the single particle relaxation rate can he simply 
found by evaluating the random potential-induced term in the action on a classical trajectory. 

As we will see now, calculation of the transport time in this formalism is much more 
elaborate. We start from the Kubo formula for the conductivity 

a = -^^{Tiv.[Gr{E) - Ga{E)]v,[Gr{E) - Ga{E)]) (14) 

where V is the system volume, is velocity operator and Gr, Ga denote retarded and 
advanced Green functions respectively. As usual, the leading contribution is given by the 
terms ~ GrGa, which can be rewritten in time representation as 

/.QO pQQ n 

a = — dTi y^ dT2 J d^R{v.GR{0, R- Ti)h,Gn(i?, 0; -T^)) (15) 

The product of the Green function can be expressed in terms of the path integral, in analogy 
with eq.(pdj) (we omit here the vertex velocity operators for simplicity and restore them in 
the end of calculation): 



(Gfi(0,/?;Ti)GA(i?,0;-T2)) 

m(Ti)=R 



/ri(0)=0 

1 pTi rTi 

2 Jo Jo 

rTi rT2 



rri[li)=n rr-i{T2)=R I r 

/ Dr:i_{t) Dr2{t)exp\i 

Jn( 0)=0 Jr2(0)=0 \ Jo 



mr'i^ 



dt- 



— I 



f'^2 mr'2^ 

dt 



dt(iFfF[ri(t) — Ti(F)] — 



1 rT2 j‘T2 

2 Jo Jo 



dt dt' W[r2it) — r2{t')] 
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dt dt' W[ri{t) — r 2 {t')] 



(16) 



One may attempt to evaluate eq.(pTD by the same simple saddle-point method, which led 
us from eq.([Tl|) to eq.(lT^). The saddle point trajectory for the action in (^) is 



5 



(17) 



n'(i) = i^R ; 
1 



r2‘{t) = 

J-2 



However, this does not give the leading contribution to (|I6D . In mathematical terms, the 
failure of the usual saddle point approximation is related to the fact that for Ti and T 2 



close to each other (this condition is implied by the factor exp{iE{Ti — T 2 )} in eq. ([T5|) ) the 
impurity-induced terms in the action in ( |T^ ) nearly cancel each other on this trajectory. 
Moreover, this happens always, when Vi and V 2 traverse the same trajectory in the same 
time, so that the vicinities of all such trajectories may be expected to be relevant. Physically 
speaking, this is because the conductivity is determined by the time scale in which the 
direction of velocity of the particle changes by an angle of order tt. Since in a long-range 
potential small-angle scattering is typical, this happens after many scattering events. The 
corresponding classical trajectory is smooth, but globally is typically very different from 
a straight line. These considerations suggest to study contributions from a class of paths 
ri(t) and r* 2 (^), which fluctuate by a small amount about a common smooth, but otherwise 
arbitrary, trajectory pit), traversed with a constant (by absolute value) velocity vp. Though 
these trajectories are not really saddle points of the action, they are “nearly saddle points” 
and can be expected to dominate the path integral. Thus, we represent the two paths T’i(t) 
and r 2 {t) in (pT|) in the form 



riih) = pit) + ; 

^ 2 (^ 2 ) = Pit) - \r-it) ; 
Tj + 1-/^, 



ti = t — = t 



tj 



T 2 tp — /2 



(18) 



where we define = (Ti + T2)/2 and t_ = T\ — T 2 . The variable t in eq.(p!8D is confined to 
the interval [0,t+]. Note that the typical values of and are expected to be ~ 1/T^ 
and ~ Ttr, so that -C t+. Now we expand the action (|^) in r’-(t). The leading 
contribution from the kinetic energy terms reads: 



Sk.in. ^ 



rt++t./2 



m 



/o 
im t- 






ft+-t-/2 



dt 



m 



^2 






= / dtp + dSkin 

I r_i_ JO 






dSkin ^ 

For the disorder-induced terms in the action 



ira I dtp(t)r_(t) 

JO 



(19) 



we use the Taylor expansion 



Wi\r + <5r|) ~ IF(r) + diWir)6ri + -didjWir)Sri6rj ; 

diWir) = W'ir)- ; 

r 

W'(r) 

did.Wir) = 






. \ ly' .ly' . 

• ’ ' + W"{r) ' ‘ 



( 20 ) 



with r = pit) — pit') and 5r = |[±r_(t) ± r_(T)]. Then the last three terms in (pT|) 
combine to 
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1 ft+ ft+ 



2 Jo Jo 



dtdt'W”[\p(t) — p{t')\]r_(t)r_(t') 



( 21 ) 



where we separated the fluctuations r_{t) into the longitudinal ri(t) and transverse r^{t) 
parts with respect to the direction of velocity p(t). As is expected on physical grounds 
and will be justified below, the relevant trajectories have nearly constant absolute value of 
velocity \p(t)\ = v. We also expect, in view of eq.(|[), the fluctuations r_{t) to be varying 



slowly on the time scale ^/v. Then eq.(Pl) reduces to 

1 rt+ 



6Sw — -Wq [ dtr-tity-tU) , 
V Jo 



( 22 ) 



with = — /o°° dr and assuming again fW'dr = 0. Finally, substituting (|T^) , (^2[) 
into eq. (pA|) , we get 



a = 



roo poo p C C C w 

dt+ J dt_ J d^Re^^*- J Dp{t) J Dr^{t) J Drl{t) p{0) p{t+ 

[ dtpJ — im [ dtp{t)r_{t) — -Wq [ dtr^{t)r^{t) 
Jo Jo V Jo 



47 t Jo j—oo 

f im /■*+ 

X exp < 

^1 2 A 



(23) 



We integrate first over the longitudinal fluctuations r_{t). This produces the ^-function 
n,i5(p(i)"). thus restricting the integral to the trajectories p{t) with constant velocity v = 
\p{t)\, as was expected. Further, the integral over gives then 6{E — mv‘^/2), fixing the 
value of V. Now we take the integral over r^{t): 

[ Dr_{t) exp{im [ dt[v{t)]'^r^{t) — -Wq [ dtr^{t)r^{t)} 

J Jo V Jo 



exp 



= exp 



mJvF 

2 

m Vp 

Aw^ 



dtv‘^{t) 



dt(j?{t) 



(24) 



where 0 is the polar angle of the velocity vector v. The action (|^) corresponds to the 
Fokker-Planck (diffusion) process for the angle Jr. 



mm) = 



m?v% 






(25) 



or, in discrete version. 



m 



2 wq 

'm?v% 



■6t 



( 26 ) 



The corresponding Fokker-Planck equation for the distribution function reads 
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( 27 ) 



7 / 7 “^ c) 

m^Vp ocp^ ot 



Thus, 

2 

a = e^N{E)-^ j dt+ j dcj) cos (j)P{(j),t+) , (28) 

where P{(j),t) satishes eq.(p7|) and the boundary condition P(0, 0) = S{(j)). To £x the nor- 
malization in eq.(^), we have exploited the condition of the particle number conservation: 



d^R 



dt_P^^-{GR{0,R;Ti)GA(R,0;-T2)) = 2 tiN{E). 



(29) 



The solution of eq. (p7[) has the form 



P{M 



oo 

exp 

m=—oo 



imcj) — 



m?v^p \ 

Wo J ’ 



We get therefore a = e^N{E)vpTtr/2, with 

1 ^ 

Ttr m^Vj, 



(30) 



(31) 



in precise agreement with eq.(||) 

We have seen therefore how the quasiclassical treatment of impurity potential reproduces 
the results of perturbation theory. In the next section we apply the method to the situation 
when the perturbation theory breaks down. 



III. MAGNETOOSCILLATIONS IN THE PRESENCE OF LONG-RANGE 

RANDOM POTENTIAL. 



We consider a 2D gas of charged particles subject to a uniform magnetic held B and a 
smooth random potential U{r) defined by eq.(|I]). We will assume the impurity scattering 
to be relatively weak, so that uJcTtr S> 1, where ujc = eB j me is the cyclotron frequency. Our 
quasiclassical treatment will be valid for a random potential with correlation length 






B 



(32) 



where Ip = {c/eByP is the magnetic length. Eq . (|3^) is opposite to the condition of ap- 

The de Haas-van Alphen 

by 



plicability of the self-consistent Born approximation 



oscillations (dHvAO) of the density of states (DOS) were studied in this regime in 
approximate summation of all orders of perturbation theory. We will demonstrate that this 
can be achieved in a more elegant way from the path integral formalism. We will also show 
that the Shubnikov-de Haas oscillations (SdHO) of conductivity can be described in this 
way, as well. 

We consider hrst the DOS, which can be found from the single-particle Green function 
(0) as 
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( 33 ) 



p(E) = -ilmGH(B) ; 

7T 

poo 

Gn{E)= {Gn{0,0-T))E^^dT 
Jo 



In the quasiclassical approximation, the Green function Gr{E) can be represented as a sum 
over closed classical orbits m 



Gr(E) 



m 

~2 



-ln(-E) 

7T 



Y, D 13 exp i 80(E) i , 

0 ] 



(34) 



where (3 labels the orbits, is a factor originating from the path integration over the 

vicinity of the classical orbit /3, and 80 (E) is an action in the energy representation. In the 
absence of the random potential, we would have just a free particle in uniform magnetic held. 
The classical trajectories are then the cyclotron orbits with radius Rc = vpjiOc- They can 
be classihed by the winding number k, specifying the number of times the orbit is traversed. 
Eq.(P4D takes then the form 



Gf(E) 



m 

~2 



1 

7T 



In(-E) - 2 iY exp 

k=l 




(35) 



where 0(E) is the step function. In particular, it is easy to check by using the Poisson 
resummation formula that (|35|) gives the correct expression for the DOS in terms of the sum 
over Landau levels: 



-| 00 

P(E)^TnTT.>>\E-^c(N+ll2)] 



(36) 



Since the impurity scattering is assumed to be weak, we neglect its inhuence on the classical 
trajectories and on the prefactor D 0 , in full analogy with calculation of in Section 2. This 
gives, instead of (|35|) , 



Gr(E) 



m 

y 



1 

71 



H-E)- 2 ,Y. exp 

k=l 



1 27iki 


\E 11 

^ 0 


- 8 wk‘^]o(E) 


{ 


-Ct^c 2. 


J 



(37) 



where 8 w is given by the second term in the action in (|Tl|) evaluated on a cyclotron orbit 
of winding number k = 1. It is easily found to be equal to 

dr j) dr'W(\r — r'\) 

= ^J(dq)W(q)[27iRMqRp)f 

dqqW(q)jQ(qRc) (38) 





If the correlation length ^ of the random potential satishes the condition ^ <C Rc, 8 \y takes 
the form 



9 



1 roo ^ 27T 7T 

Sw = / dqW(q) = / drW(r) = , 

UJc^F do VpUJc do OJcTs 

with Ts as found in section 1, see eq.(||). We get then for the DOS ai E > Q 



( 39 ) 



p{E) = — 
’ 2n 



/ JTJ \ / 

1 + 2 ^(-1)^ cos [2vrA;— j exp 



fc=i 



7T 



^cTfi 



( 40 ) 



or after resumming by the Poisson formula, 

OO 



p{E) = 



E 



a'‘=o ^ 



exp 



UJn 



[E-u,{N + l/2)f 



(41) 



This formula shows that the Landau levels acquire a Gaussian form. They are well resolved 
if uJcTs 3> 1. In the opposite case, uJcTs "C 1, the representation (|40|) is appropriate, where 
all harmonics except the hrst one can be omitted: 



p{E) =m/27i + Posc{E) ; 



pose - cos{2ttE / ujc)e 

7T 



(42) 



The dependence of the amplitude of oscillations on the magnetic field has the same form 
~ exp(— vr/cUcTs) as for the short range potential, so that if one plots log posc versus 1/B 
(so-called Dingle plot |^), one expects to get a linear behavior. 

In the case of ultra-long-range potential with ^ S> i?c, we hnd Sw = 27r^lT(0)/a;^. The 
Landau levels have again a Gaussian shape: 



P{E) = 



at=o 



^2-kW (0) exp 



2W(0) 



[E - o^,(iV + 1/2)]' 



(43) 



which in this case reflects their inhomogeneous broadening. The oscillating part of the DOS 



in eq.(h^, Posa is now equal to 



~ --cos(27TE/a;,)e-2"''^(°)/‘^' , 



7T 



(44) 



so that the Dingle plot is expected to be quadratic. 

As we have already mentioned, the above results for the DOS oscillations in a long 
range potential were obtained in [Q by resummation of the perturbation theory expansion 
0. Besides being simpler and physically more transparent, the present derivation has the 
advantage that it can be straightforwardly generalized to SdHO of conductivity or to the 
case of random magnetic held. 

The conductivity oscillations were recently considered in the framework of the path 
integral approach for the case of short-range random potential (r^^ = t^) in p3|,|2^. The 
authors of these papers started from the Kubo formula (pA|), representing each Green function 
as a sum over classical trajectories in uniform magnetic held and taking scattering into 
account by including the factor exp(— t/r^). A trajectory is characterized by a number of 
cyclotron revolutions k. Therefore, eq.([l5|) is reduced to a double sum over the winding 



numbers fc/j, k^- The non-oscillating contribution to the conductivity corresponds then to 
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the terms with kji = kA, whereas the n-th harmonic of the oscillations is described by the 
terms with = n. This holds for our problem of smooth random potential as well. 

The only difference is that the relevant trajectories are now not exactly cyclotron circles, 
but rather are very close to them, however drifting a little from one revolution to another. 
The typical total number of revolutions of each trajectory is of order of u>cTtr 3> 1; in this 
time the shift of the center of cyclotron movement for a typical trajectory is of order Rc. 
This is completely analogous to the consideration of conductivity in zero magnetic held in 
section 2, where the characteristic trajectories dominating the path integral were smooth 
(i.e. locally close to straight lines), but not really straight lines. 

The non-oscillating contribution (Tq to the conductivity is given by the trajectories for 
Gr and Ga having equal number of cyclotron revolutions and following closely each other 
within a distance -C as in section 2. Then the impurity-induced terms in (|^) cancel each 
other again to a great extent, leading to Ttr 3> r^. The exact evaluation of the conductivity 
(To is however much more complicated in the present situation, since the trajectories return 
to positions close to the preceding ones after 1 , 2 ,... revolutions and may interact through 
the impurity correlator. This may lead to a deviation of Ttr from its value calculated in 
zero magnetic held in section 2. We leave this problem aside in the present article and 
concentrate on the exponential damping of the oscillations by disorder. 

The leading contribution to oscillations (the hrst harmonic) corresponds to the case when 
one of the trajectories has one extra cyclotron revolution as compared to another one. Then 
the contribution to the effective action in (^) from this part of the trajectory is not cancelled 
and leads to a suppression of the amplitude of oscillations. We get 

a = ao + aiCos{2TE /ujc)e~^^ + ■ ■ ■ , (45) 



with Sw given by eq. (|38|) . Thus, impurity scattering leads to the same e exponential 
damping of the SdHO, as for the DOS oscillations. Depending on the relation between ^ 



and this factor takes the form 



^-Sw 



^ — TTjljJrTs 
,-27r2W(0)/o;2 



Rc 



(46) 



as in eqs. (^, (^. 



IV. CONDUCTIVITY AND MAGNETOOSCILLATIONS IN RANDOM 

MAGNETIC FIELD. 



In this section we consider the oscillations of DOS and conductivity in the case when 
the carriers are subject to a random magnetic held rather than to a random potential. The 
properties of a quantum particle moving in a random magnetic held have been intensively 
discussed in the recent literature ||T^ , ^5| -^. We are aware by now of three completely dif- 
ferent physical systems, for which the random magnetic held problem is relevant. In the 
hrst class of system the random magnetic held is generated by a substrate on top of which 
the 2D electron gas is placed. Realizations of the substrate are a type II superconductor 
with randomly pinned hux lines [|l|, or a demagnetized ferromagnet with randomly oriented 
domains [0]. The second is the state with spin-charge separation of high-T^ superconducting 
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materials p-§[. To describe it, one introduces a fictitious U{1) gauge field interacting with 
charge carriers. Here, the transverse (magnetic) component of the gauge field is the most 
important. The existence of gauge field fluctuations in high-T^ compounds can be inferred 
0 from experimental observations of unusual weak localization corrections to the magneto- 
conductance 0. Finally, the third application concerns the FQHE system in the vicinity of 
z/ = 1/2 filling factor of the Landau level. This will be studied in more detail in the next 
section. 

We consider a Gaussian distributed random magnetic field h{r) (perpendicular to the 
2D plane) with the correlator 



{h{r)h{r')) = T6^'^\r — r') 



(47) 



We will assume the weak disorder case, which means F -C mEp. Let us first reconsider 
the evaluation of scattering rates in section 2 for this type of disorder. For the one-particle 
Green function we have, instead of (|np: 



(G«(0,i?;T)) 





dt 



mr'^ 

~Y~ 
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dt dt' r j {t)rj (t'){ai{r (t) )aj{r{t'))) 



(48) 



where a,{r) is a vector potential corresponding to the magnetic field h{r). However, the 
correlator (a^aj) is gauge-dependent, so that eg. (|48|) does not allow to define meaningfully 
the single particle relaxation time r^, as it was possible for the random potential case. We 
will discuss the problem of definiton of below. 

For the transport time, the problem of gauge invariance does not apply. The analogue 
of eg. (0) reads 



{Gr{0,R;T,)Ga{R,0;-T2)) 



rrx{Ti)=R rr2{T2)=R f pTi rnV-i'^ f 

= Dri{t) Dr 2 {t)exp\i dt— i 

7ri(0)=0 Jr2(0)=0 ( 4o 2 Jo 2 

Jo Jq dtdt'Gi{t)rij{t'){ai{ri{t))aj{ri{t'))) 

~2\^) lo lo ^^^^'^2i(t)f'2i(^')(«i(’’2(^))ai(^’2(^'))) 

+ (^) ^ ^ dtdt' ru{t)r2j{t'){ai{ri{t))aj{r2{t')))\ 

ri{Ti)=R ^ pr2(T2)=R _ _ , f pTi TflVi^ f'^'^ TTIV^ 



mri'’' /■^2 rnr 
— i dt 



^2 



Dri{t) I Dr 2 {t) exp [ dt 

/ri(0)=0 Jr2{0)=0 [ JO 



2 -'L 



where we used the Stokes theorem to rewrite the effective action in the explicitly gauge- 
invariant form. Here Sno is so-called non-oriented, or Amperean, area inside the closed curve 
formed by the two trajectories Ti(t) and r 2 {t) We recall, in particular, that the non- 



oriented area enclosed by k windings of a trajectory is the geometrical area multiplied by 
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k'^. We repeat now the derivation of Ttr performed in section 2. For the disorder- induced 
part of the action we hnd, instead of (^2[), 



/■*+ 

= ^ dt\r^{t)\ 

2 Jo 



(50) 



Integral over the fluctuations eq.(|2^), now has the following form: 

J Dr^(t) exp{—im j dt[i)\'^r^ dt\r_\} 



~ n y exp 

Sti 



n 



Fn 6tj 



5u [0^v/2)6ti]‘^ + m^[6v{ti)]'^ 

This implies that the scattering angle ~ [5v\-^ /v obeys the Cauchy distribution 

PU)S] = 1 (F/2m)dt 

7T(F/2m)2(dt)2 + (d0)2 

The Boltzmann equation corresponding to eq.(|52D is found to be 

dP{4>, t) 



(51) 



(52) 



m 



= J d<l>'w{<p - 4>')IP{4>', t) - P(>l>, i)l , 



(53) 



with 



1 F 1 

^ (</>) = -Ti—-! 1 0<vr (54) 

7T 2m (p^ 

This can be easily checked by solving eg. (|53|) by means of the Fourier transform in the 0- 
space. The quasiclassical method describes correctly only the small angle scattering, so it 
is able to give an expression for the differential scattering rate w{(p) for <p ir only. The 
transport scattering rate l/rtr = — cos cj))d(f) can be found in this way up to a 

numerical coefficient only: l/rf^ ~ F/m. The exact value can be found from the perturbation 
theory calculation, which gives [^] 



W(0) 



F 

Svrm 



cot^ <p/2 , 



in full agreement with eq. (|54|) . Consequently, 



1/Ttr = F/4m 



(55) 



(56) 



As we can see from eqs. (|5^) , (^), the total scattering rate l/r* = / dcpw^cp) diverges at 
0 — >■ 0. This is related to the fact that the contribution Sh to the action from the random 
magnetic held is proportional to the area, rather than to the length, of the trajectory. 
We will see below, when studying dHvAO and SdHO, that the cyclotron motion of the 
particle provides a natural regularization of this divergency. However, this regularization 
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is determined by the geometry of the experiment considered. Thus, the single particle 
relaxation time in the random magnetic field is dependent on the geometry of the problem. 

We turn now to the consideration of magnetooscillations in the case when a uniform 
magnetic held B is applied in addition to the random one. As in section 3, we will assume a 
relatively strong held, meaning oJcTtr 3> 1, or according to eg. (|56D , mojc S> T. In full analogy 
with eq.(^), the DOS can then be written as 



P(E) 



m 



°° / /? \ 

1 + 2 ^ (—1)^ cos ( 27rk — ) exp (—k'^Sh) 
k=i ^ W/ 



(57) 



with Sh being the random magnetic held action on a simple (winding number 1) cyclotron 
orbit. For a 5-hke correlated magnetic held, eq.(^7|), it is equal to 

Sh = ^Tsno=^'TTrRl = irET/mu^ (58) 

Resumming eq. (|^ with the help of the Poisson formula, we hnd 



P{E) 



1 

27llfi 







(59) 



As in eg. (^I|) , the Landau levels have a Gaussian shape. However, in contrast to the case 
of random potential, their width does not increase with the magnetic held, but is instead 
proportional to When the oscillations are relatively weak, they can be characterized 

by the hrst harmonic with an amplitude 



pose 






m 

cos 

71 







(60) 



These results for the amplitude of oscillations can be generalized to the SdHO of conduc- 
tivity, as in the preceding section. We get again eq. (|45|) , with Sw replaced by Sh, eq.(^). 
We can also consider a random magnetic held with hnite correlation length ^ and correlator 



{h{r)h{r')) = U{\r — r'\) 



(61) 



We hnd then 



Pose, OSC 0 ^ 6 ) 



Sh = 7^ d^r dVH(lr-r'l) 

2 J\r\<Rc J\r'\<Rc 



= 7iRl r "^U{q)Jl{qR, 
0 



7rEJd^rU{r) 

mul 

27t‘^E‘^U{R) 



Rc<^f, 



(62) 



where U{q) is the Fourier transform of U{r). Therefore, the Dingle plot will be quadratic 
for ^ and quartic in the opposite case. 
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V. MAGNETOOSCILLATIONS NEAR THE i/ = 1/2 FILLING FACTOR OF THE 

LANDAU LEVEL. 



In this section, we consider a realistic model describing the electron gas in GaAs — 
—AlGaAs heterostructnres where the FQHE is observed. The system is formed by the 
2D electron gas of density rie and by the positively charged imparities located in a layer 
separated by a large distance ds from the electron plane. Each imparity creates a potential 
of the form 

J ; Vo{q) = (63) 

where Vi is the projection of the impurity position to the 2D plane and e is the dielectric 
constant. We briefly consider first the magnetooscillations in this system in low magnetic 
fields. The charge carriers are then characterized by their lattice mass mb, which for GaAs 
is equal to rrih ~ 0.07me, mg being the free electron mass. The potential Vo{q), eq.(|6^), gets 
screened by the 2D electron gas into 



v{q) = — 
rrib 



(64) 



When writing eq . (p4f) , we assumed that the DOS determining the screening is practically 
constant: p ~ mb/27r. Taking into account the magnetooscillations of the DOS here would 
lead to a non-linear screening, see 



37 



However, since our primary interest is in the region 
of large p where the amplitude of oscillations is small, we can neglect this non-linear effect. 
We also assume the impurity positions to be uncorrelated (see, however, the discussion in the 
next section). Then we hnd that the total random potential of the impurities is described 
by the correlator IU(r), see eq.(|l|), with the Fourier transform 



W{q) = Hi f— ^ e 
KrUbJ 



(65) 



According to section 3, we hnd then for the amplitude of oscillations 



cos 



"27T^neC' 

eB . 



■.-Sw 



(66) 



with Sw given by eq. (|38|) . In particular, for ^ dg we hnd 



pose, oc - cos 



'27r^neC^ 



eB 



To TTVf 



dqW{q) = 



Hi 



exp(-7r/o;gro) ; 



27ry/2 



TUbdg \rieJ 



whereas for Rr -C d. 



Pose, (^xx - COS 

1 
T-* 

A) 



( 2TT^neC\ 




eB 


- (ttMto) 



Tflhds 



(67) 



( 68 ) 
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For the systems under consideration the usual assumption is that concentrations of donors 
and charge carriers coincide: rze = Uj. Then the condition of weak oscillations oOc^o -C 1 
reduces to Rc 3> dg- Therefore, in the region of small oscillations (where our derivation is 
justihed), eq.(|67|) has to be applied. 

Now we consider the oscillations in strong magnetic held, near half hlling of the Landau 
level: u = 27rcne/eB ~ 1/2. It was observed experimentally that in this region the longitu- 
dinal resisitivity shows oscillations with minima ai u = p/{2p ± 1), very much reminiscent 
to its behavior in low magnetic helds where conventional SdHO take place. This feature did 
not hnd an explanation within the original hierarchy theory of the FQHE P5| , PU| . To explain 
it, Jain Q] proposed a concept based on converting the electrons into composite fermions 
by attaching to them an even number of hux quanta. A held-theoretical formalism based 
on the Jain’s idea was developed by Lopez and Fradkin 0. In this approach, the statis- 
tical transformation of electrons into composite fermions is implemented by introducing a 
Chern-Simons (CS) gauge held interacting with electrons. 

Following a similar approach Halperin, Lee and Read [jT^ developed a theory for the half 
hlled Landau level (see also |^). This theory gives an explanation for many experimentally 
observed properties of the u = 1/2 state, such as a non-zero value of the longitudinal 

M, and a dimensional 



resistivity, an anomaly in the surface acoustic wave propagation 
resonance of the composite fermions [|1]. It predicts the formation, at half hlling, of a 
metallic state with well dehned Fermi surface. From this point of view, the v = pj {2p ± 1) 
series can be considered as the usual u = p Shubnikov-de Haas oscillations (SDHO) for the 
composite fermions, providing an explanation for the prominence of the above FQHE states. 

Details of the CS gauge held formalism can be found e.g. in P,pR[| and are not presented 
here. The statistical transformation attaches to each electron an even number 0 of hux 
quanta of the CS gauge held. To describe the vicinity of the z/ = 1/2 state, we take 0 = 2; 
the same formalism with 0 = 4 can be applied to the = 1/4 state. In the mean held 
approximation, the statistical magnetic held B 1/2 = A.-Kcrteje cancels exactly the externally 
applied held B at u = 1/2. When the hlling factor u is tuned away from u = 1/2, the 
ehective uniform magnetic held is equal to Re// = B — R 1 / 2 . For u close to 1/2, the number 
of hlled Landau levels of composite fermions p 3> 1, so that the problem can be considered 
quasiclassically. 

Although a static impurity creates a scalar potential ( |63D only, it acquires also a vector 
component due to screening by fermions and mixing with the CS held. In the random phase 
arrroximation one gets 



= (69) 

where we united scalar oq and vector a potentials in a covariant vector a^; the vector 
represents the bare impurity potential, eq.(pB|) and therefore has only p = 0 non- zero 
component. The tensors and represent the bare gauge held propagator and the 
current- density response tensor of the composite fermions, respectively. 

To evaluate eq.(^) we use the Coulomb gauge diva = 0, go to the momentum space and 
choose the momentum q to be directed along the x-axis: qx = q, Qy = 0. Then has only 2 
non-zero components corresponding to /i = 0, y, and both K and U become 2x2 matrices 
0 : 
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U,M) = 
= 



( -171*1271 -iqCTxy 

\ iq<Jxy XQ^ - ‘^i^rie/qkp 

f v{q) 27ii^/q\ 

y —27Ti(f)/q 0 j 

/ Vo{q)e~^^^i \ 



(70) 



where m* is the effective mass of fermions, x = l/247rm* is the magnetic snsceptibility, 
v{q) = 27re^/(eg) is the Conlomb propagator, and axy is the Hall condnctivity of composite 
fermions. 

Snbstitnting (|70|) in (§^), we find 



a^^{Q) = 



Vo{q)e 

+ (0s + 1)2 + g V ^M*/q 



(71) 



where s = 27raxy ~ p in the limit uicTt S> 1 . 

Depending on relations between parameters of the problem, one can find varions regimes 
of behavior of the oscillations amplitnde. We will concentrate on a regime which is the most 
relevant to the experiment. Let us compare the first and the second term in denominator of 
(0)- As we will see below, the typical momenta are q ~ (2ds) and we get for 0 = 2 



m*v{q)/27i 



^ 2 
m e 

Aeqs'^ 



m*e^ kpds 
ekp 2 p 2 p 



50 
~2 ’ 



(72) 



where kp = A/4vm^, and we used typical experimental parameters Ug = 1.1 ■ 10^^cm“2^ 
dg = 80nm, and the experimentally estimated value for the ratio m*e^ / {ekp) ~ 10. For the 
not too large p we are interested in, it is thus a reasonable approximation to neglect all but 
the first term in the denominator of (|7l|) . This gives 






27T 

m* 




-qda 




(73) 



The random field action Sr (analogous to Sw in section 3 or Sh in section 4) is given by 

5'r = ^ ( (^ ^ = aodt - j adr^ ^ , (74) 

where the integration goes around a cyclotron orbit. Averaging over the impurity configu- 
rations, we find 



Sr = (27T0)^nj 



X 



^-2qds 



kp 



= nj(4vrV-Rc)" 



^ 2,qds 



-^Jo{qRc) H — JiidRc 

kp q 



(75) 
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Here f dl means integration along the cyclotron orbit and corresponds to the electric held 
contribution, whereas / d?r goes over the area surrounded by the orbit and describes the 
magnetic held contribution. Taking into account that we have Rc/2dg = 

p/^Aimedl p/10 ^ 1. Thus for relevant momenta q ~ 1/(2^^) and level numbers p, 
qRc -C 1 is a reasonable approximation. In this case eq.([75D reduces to 



Sr 




Hi / 27me 

Ue Uedl \m*UJc 



4 



(76) 



Note that electric and magnetic held huctuations give equal contributions in this limit. 
According to sections 3,4, this gives for the oscillating part of the conductivity: 






cos 






eB, 



eff 



exp 




(77) 



where we introduced a parameter r /^2 which is given according to eq.( ff6|) by 



' 1/2 — 



m 

T 



dl 



1/4 






(78) 



Let us briehy consider now ehect of hnite temperature T. First of all, the SdHO are then 
suppressed by the usual factor Dt = /ujc) j sinh(27r^T/a;c) originating from the Fermi 

distribution |^. In addition, the fermions are scattered by the thermal huctuations of the 
gauge held. The propagator of gauge held huctuations is given by 






-1 



(79) 



In particular for the Du component determining the magnetic held huctuations, we get 



Dii{q,u;) = {2iume/qkF - XQ 



2\-l. 



X = 



2Tim* 



1 


( 


— — + 
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V 07 J 



+ 



v{q) 

(27T0)" 



In the quasistatic approximation we hnd 



f duj2T T 

{AiAi)q = / ImHii = — 

J 2n oj XT 

and consequently for the amplitude of magnetic held huctuations 



(80) 



( 81 ) 



(A/!), = T/x. 



(82) 



If uJcTt S> 1, we have s ~ p and y = 12p^y. Therefore, the ehective magnetic susceptibility 
X much exceeds its bare value y, leading to a strong suppression of the huctuations ( 0 . 
(^). A similar suppression of the gauge held huctuations in external magnetic held was 
found in [H where the magnetoconductivity of the doped Mott insulators was studied. The 



contribution of the huctuations (|8^) to the random held action 5^, eq.(^), is 



18 



27t2 T 



( 83 ) 



~ , 

X P 

i.e. is small at p S> 1 compared to the standard term — \\\Dt and decreases with p. 

The slope of the dependence of Inpo^c on 1/5 at T S> o;c is conventionally used to 
extract a value of the effective mass m* [^. This procedure is based on the assumption 
that Pose oc Dt, so that Inpo^c = 27i^Tm*c/eB + const. The additional attenuation factor 
exp(— eq.(^), would then lead to a hetitious 1/p correction to the effective mass 
resulting in its apparent decrease with p at moderately large p. 

Let us compare our hndings with available experimental results. In Fig.l we present low- 
temperature experimental data for the amplitude of p°^‘^ from (T = 0.19iL, Beff >0). It 
is seen that they can be htted well by exp [—(pK/ujcT*Y], whereas a conventional exp(— 7r/o;cT) 
£t is much worse. Therefore the data apparently show the behavior In po^c oc l/o;/ predicted 
by eq.(^). The value of the parameter which is found from such a £t is r /^2 = 16- 10“^^s 
At the same time the theoretical estimate according to eq. ([78|) (with use of the parameters 
of [^) gives r/y 2 — 2.4 ■ 10“^^s if one uses the experimental value of the effective mass 
m* = 0.7me. A similar discrepancy is found for the low-held relaxation time: eq.(|67|) for 



mb = 0.07me gives rj ~ 0.6 ■ 10“^^s, whereas the value quoted in is tq — 9 ■ 10“^^s. We 
note also that the theoretically estimated values for the transport relaxation rate at i/ = 1/2 
are typically 4 times greater than extracted from experimental mobilities Therefore 



the theory seems to overestimate relaxation rates systematically. This situation has been 
discussed previously The considerable increase of relaxation times was attributed to 



the correlations in positions of charged impurities due to their mutual Coulomb interaction 
1^4611 . We still encounter a problem, however, when trying to explain the above discrepancy 
in the value of in this way. The correlations between impurities lead to a damping 
of the correlator W{q), eq.(|65|) by a certain g-dependent factor. Since the mechanism of 
suppression of SdHO and the characteristic momenta are the same near 5 = 0 and z/ = 1/2, 
we expect a suppression of the action Sw in (|66D and Sr in (0) by roughly the same factor. 
For the case of low helds this factor is To{exp.) /To{theor.) ~ 15. At the same time, to 
reconcile the experimental data in the vicinity of z/ = 1/2 with eq.(]78D, we need this factor 
to be [r// 2 (ea;p.)/r //2 (theor.)]'^ ~ 2300, i.e 150 times larger! It is not clear to us, what could 
be the source of such a drastic weakening of the random helds. 

More recently, experimental data on SdHO near u = 1/2 on a better quality sample 
have been published. The obtained Dingle plot [Fig. 3(a) of Ref. [^] is again highly nonlinear 
and can be well htted by our formula ([77|). The ht yields in this case the value r /^2 = 
9.5 ■ 10“^^s, whereas the theoretical estimate according to eq.(^) using the appropriate 



values of Ue,ds and m* gives r /^2 = 2.0 ■ 10“^^s. The discrepancy is somewhat smaller than 
for the sample from but still very large. 

Now we discuss experimental data at higher temperatures. As a result of their analysis, 
an experimental value of the ehective mass as a function of p was obtained in |]4^ , |47| , ^ . 
At moderately large p m* was found to be slowly decreasing with p in agreement with our 



results. For larger p, a sharp increase of m* was observed in |^8||, the origin of which is not 
clear to us. Let us note that the suggestion made in to explain the non-linearity of the 
Dingle plot by this variation of the ehective mass is not supported by our results, since m* 
drops out from eq. dZ3)- 



19 



VI. CONCLUSIONS. 



We have shown in this article that the path integral formalism in the quasiclassical 
approximation allows one to study the density of states and the conductivity of a system of 
charge carriers scattered by a long-range correlated random potential or a random magnetic 
held. We focussed attention on the oscillatory behavior of these quantities as a function 
of an applied magnetic held. For a random potential of correlation length we found the 
magnetooscillations in both quantities to be attenuated exponentially, the exponent being 
oc for Rc ^ ^ and oc for Rc where Uc is the cyclotron frequency and Rc the 
cyclotron radius in the external magnetic held. In the case of a random magnetic held of 
correlation length ^b, the magnetooscillations are again exponentially damped, this time 
with exponent oc for Rc 3> and oc for Rc -C ^b- Our results show that the 
amplitude of magnetooscillations as a function of the magnetic held can be used to identify 
the scattering potential or random magnetic held. 

Finally, we considered the FQHE system near z/ = 1/2 in the composite fermion picture. 
Here the Chern-Simons held of the hux tubes attached to every fermion gives rise to static 
random held huctuations at impurity sites, where fermions may be trapped. For experi- 
mentally relevant values of input parameters, the damping of magnetooscillations is found 
to be proportional to exp[— (7 t/ci;cT// 2)"^]- This is in good agreement with experimentally ob- 
tained highly non-linear Dingle plots which can be very well htted by a quartic dependence. 
However, the experimental value of the parameter is about 8 times greater than our 
theoretical estimate. This means that the random helds are much weaker than we expect 
them to be. Taking into account correlations between impurities seems not to resolve this 
discrepancy, as we discussed in the end of the preceding section. It is not clear to us at 
present whether this additional weakening of random helds can be explained within the 
composite fermions theory or else implies a certain inconsistency of this theory. 

This work was initiated by the late Arkady Aronov, and was partially done in collabora- 
tion with him. We very much regret that he did not live to see the completion of the work. 
We are grateful to Yehoshua Levinson for useful discussions. This work was supported by 
the Alexander von Humboldt Stiftung (A.D.M.), SFB 195 der Deutschen Forschungsgemein- 
schaft (A.D.M. and P.W.) and by the German-Israel Foundation for Research (E.A.) 
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FIGURES 



FIG. 1. Dingle plot. Logarithm of normalized amplitude of resistivity oscillations 
/ ^p) , with Dt = sinh(27r^r/wc)/(2vr^r/u;c), as a function of inverse effective magnetic 



held Rg//- Experimental data from 
hts with eq.(77). 



(squares) and (circles) are presented, as well as their 
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